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Abstract 

We report computations of the short-distance and the long-distance (scahng) contributions to 
the square-lattice Ising susceptibility in zero field close to Tc. Both computations rely on the use 
of nonlinear partial difference equations for the correlation functions. By summing the correlation 
functions, we give an algorithm of complexity 0(A^^) for the determination of the first A^ series 
coefficients. Consequently, we have generated and analysed series of length several hundred terms, 
generated in about 100 hours on an obsolete workstation. In terms of a temperature variable, r, 
linear in T/T^ — 1, the short-distance terms are shown to have the form rP(ln|r|)'^ with p > q^. 
To 0(r^^) the long-distance part divided by the leading r~^'^ singularity contains only integer 
powers of r. The presence of irrelevant variables in the scaling function is clearly evident, with 
contributions of distinct character at leading orders Irp^ and |t|^'^'^ being identified. 
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I. INTRODUCTION 

The two-dimensional Ising model has been extremely useful as a testing ground for new 
theoretical ideas and methods in the study of phase transitions and critical phenomena. Our 
present understanding is the result of a series of dramatic developments spanning more than 
half a century, starting with Onsager's exact computation of the free energy ]l|, followed by 
Yang's derivation of the spontaneous magnetization [Q] and by the work of many researchers 
on the correlation functions, including Toeplitz determinantal formulae [Q, exact expressions 
for their behaviour at large separation [^ , and nonlinear partial difference equations for their 
efficient computation [§, ^, 0], to mention only those results which are used in the present 
work. All results above apply to the zero-field case. While an exact expression for the 
susceptibility as the sum of two-point correlation functions over all separations [Q exists, a 
useful closed form expression does not. Moreover, as we discuss, there are strong indications 
that the susceptibility has a natural boundary in the complex plane §, ^, a feature which 
rules out any expression in terms of the "standard" functions of mathematical physics. 

Nevertheless, it is desirable to obtain as detailed information about the susceptibility as 
possible, not only because of its physical importance, but also because of the significant 
role it plays in ideas about scaling and the renormalization group. In the vicinity of the 
ferromagnetic critical point at temperature T = T^, the susceptibility exhibits a singularity 
of the form 

p-\± = Co^{2K,V2Y^'\T\-'/'F^{r) + B{r). (1) 

Here /? = {ksT)-^, r = \{s-^ - s), s = sinh2i^ and 2Kc = ln(l + ^2) with K = (3J 
the conventional Ising model coupling constant. The scaling-amplitude functions F±{t) are 
normalized to unity at r = 0. As a consequence of the exact knowledge of the long-range 
correlations, the coefficients Cq^ were calculated exactly 0] in terms of the solution of a 
Painleve III equation. Additionally, the leading behaviour of both F±{t) was computed 
to be 1 -|- ir. The antiferromagnetic susceptibility, on the other hand, is dominated by 
the short-distance correlation functions and has leading singularity (const x rln|r|). Such 
short- distance "background" terms are present as well in the ferromagnetic susceptibility 
and are denoted by B{t) in (|l|). The leading amplitudes of the analytic and singular parts 



of B{t) were computed for a general wavevector dependent susceptibility in [O, O] 



An analysis [|T3| of a 51 term high-temperature series by means of differential approxi- 



mants yielded two further correction terms in the scaling- amplitude function F+, with nu- 
merical amplitudes close to rational values, |r^ + jqT^, and confirmed that the same scaling- 
amplitude function is numerically consistent with the first 11 terms in the low-temperature 
expansion. These results agreed with the prediction |T^ that the corrections to scaling are 



entirely due to the nonlinearity of the scaling fields and not to the presence of irrelevant 
operators [Q. However, a recent analysis of 115 term high- and low-temperature series ||T5[ 



showed that this prediction appears to break down in the amplitude of r^. 

The study reported in this letter substantially improves on all the above results. We 



extend the methods of |]TI], |1^] to compute both antiferromagnetic and ferromagnetic back- 
ground amplitudes on the isotropic lattice to 0(r^^(ln |t|)^). All such terms are seen to be 
of the form rP(ln Irl)"^ with p > q^. We simultaneously compute high-temperature series to 
order 323 and low-temperature series to order 646 in 123 hours on a 500MHz DEC Alpha 
with 21164 processor running Maple"""^ V version 5.1. 

We analyze these series by two independent methods, making use of the computed back- 



ground amplitudes and the known complex singularity structure ||^, 15 to obtain the scaling- 
amplitude functions F± to 0(r^^). 

Several important conclusions can be drawn from our results. Firstly, only pure integer 
powers of r enter the scaling- amplitude functions and no logarithmic terms are present. 
Secondly, the high- and low-temperature scaling-amplitude functions are not equal to each 
other. The amplitudes start to differ at O(r^). Thirdly, the amplitudes of r"^ and r^, which 



are clearly rational, are not those predicted by simple two- variable scaling [T^. We surmise 



that at least two irrelevant operators must be invoked to account for the above results — one 
entering at r^, the other at r^. 

Further remarks on the scaling implications of our work can be found in section while 
the remainder of the letter will outline the methods by which the ferromagnetic results were 
obtained. A fuller account, including details of the antiferromagnetic singularity, will appear 
elsewhere. 



II. SINGULARITY STRUCTURE AND NATURAL BOUNDARY 

It was argued in 0] that on the anisotropic lattice, the contribution to the susceptibihty 
of the high-temperature graphs with 2N vertical bonds contains more and more poles as A^ 
increases, and that in the limit N —^ oo these poles form a dense set in the complex plane. 
In it was shown that in the expansion of the susceptibility in j-particle contributions 

y/-^j even "X. -L < I c, 

the higher-particle components give rise to an ever increasing number of singularities that 
appear to form a dense set on the circle |s| = 1. In fact, the two phenomena are precisely 
correlated, with the former being the highly anisotropic limit, and the latter the isotropic 
limit, of the set of singularities for the generic anisotropic model. These occur at 

cosh(2ir) cosh(2ir') - sinh(2K) cos 

J 

- sinh(2ir') cos — = (3) 

with m, m' = 1, 2, . . .j, and K, K' = pj^, pjy. It will be noted that the left-hand side of (H) 
is the denominator in the Onsager integral for the free-energy and thus we find the (to us) 
surprising result that the singularity of x^''\ ^ property of the Ising model in a magnetic 
field, is intimately connected with a property in zero field. Barring unexpected cancellation 
(and we have evidence against this) in the N -^ oo limit, we believe that this set forms a 
natural boundary. 

The existence of a natural boundary has strong implications for our series analysis. In 
the r plane, the boundary is formed on the imaginary axis by the logarithmic branch cuts 
coming from the singularities of 0, |T^. Summing the contributions of the discontinuities 
across these cuts gives Disc(/3~^x) ~ exp(— 39.76/T^) to leading order, where t = iT is 
a point on the cut. Assuming that the contribution of the physical singularity at r = 
is additive, we expect that the coefficient of t^ in the limit p ^ oo in the r-expansion 
of the susceptibility will grow as T{p/2)/aP/'^ where a = 39.76. This effect only becomes 
appreciable at order p ~ 2a and is too small to be seen in the terms we have. However it 
has the implication that the r-expansion of x can be at best asymptotic and not analytic. 
Whether this non-analyticity comes from the long-distance or the short-distance part, or a 
combination of the two, has not been determined. 



We note that the existence of a dense set of singularities also implies that the suscep- 
tibility cannot be a member of the class of functions called "D-finite, that is, solutions to 
linear differential equations with polynomial coefficients. Nevertheless, it can be shown 
from the integral expressions for the x^"''' that they individually are "D-finite (also called 
holonomic) |T^ . This phenomenon is related to that of perturbative expansions in quantum 
field theory, where any individual Feynman diagram is holonomic whereas the summed series 
may not be ||T7[| . 



III. COMPUTATION OF SHORT-DISTANCE AMPLITUDES AND HIGH- AND 
LOW-TEMPERATURE SERIES 

The essential tool for the computation of both the background amplitudes and the high- 
and low-temperature series coefficients is the set of nonlinear partial difference equations 
for the two-point correlation functions, C{m,n) = (o"ooO"m„), given in 0. These completely 
determine all the off-diagonal two-point functions once the diagonal ones {m = n) are 
given. The latter can be computed either by means of an independent set of difference 
equations or, as we have done here, directly from the Toeplitz determinant expressions. 
The susceptibility, P~'^x = XlC'(m, n) — (ctoo)^, is computed by successively adding the 
contributions of pairs of square shells Cn = Yl^i''^'^) with |?7i| + |n| = 2N and |?7i| + |n| = 
2N + 1. 

The implementation of the difference equations to obtain high- and low-temperature 
expansions is straightforward using the multiple precision integer arithmetic capabilities of 
Maple"""^ or Mathematica"""^, and the time complexity is no worse than 0(A^^). 

The key to computing the short-distance background amplitudes is to obtain expansions 
of the partial sums Sn = Yln=o ^n in t directly and to identify which terms in the series 
contribute to the short- distance part and which to the long-distance part. A combination of 
analytic work and numerical fitting leads us to a conjecture for the short- distance expansion 
of the shell sums, namely 

oo 

^sCn = N^'^ 5^(ln \Nt\Y{NtY"A%^ (4) 

p=0 

where the Ap are Taylor series in r with coefficients that are asymptotic Laurent series in 



A^ ^; the highest power of A^ multiplying r'^ in A^ is A^''. The partial sums S^ are 



N 



n=0 P=q^ 9=0 

with -R]^''' functions of A^ only. Asymptotically, for large A^, -RJ^'"^ is a sum of powers 
JY7/4+P ^ with possible multiplicative In(A^) corrections, plus a constant b^''^^ which arises 
from the small n terms in the sum (J^) where the asymptotic expressions are not valid and 
sum and integral are not synonymous. The p' are integers p' < p. 

We must assume that (^) remains valid up to A^ of the order 1/r where it can, in principle, 
be matched term by term to a large distance expansion that properly describes the roughly 
exponential exp(— A^r) decay of correlations as A^ ^ oo. Explicit matching formed the basis 



of the previous calculations of terms in the short- distance x (cf- ini, 0) but this becomes 
extremely cumbersome at higher order. Here we argue that the exponential decay implies a 
cutoff on A^ that is proportional to 1/r and that we can identify the temperature behaviour 
of terms in Sm in eqn. (j^) by the replacement A^ — ^ 1/r. All terms whose variation is as 
a fractional power of r, with possibly logarithmic multipliers, are discarded as assumed 
contributions to the long-distance part of x- Clearly all that remains is the constant part of 
i?]^''^ , namely b^''^\ and this is extracted by numerical fitting to give 

y/TsB = J2Y1 b^^'^'^^^i^^ \^\y (6) 

p=g2 q=0 

for the short-distance part of x iii sqn. ([I|). The coefficients b^P'''^ must be determined to 
very high accuracy to be useful for the subtraction procedure described in the next section; 
the complete list for p < 15 will be given elsewhere. The result p > q"^ we call the fermionic 
constraint since it can be traced back to the Toeplitz determinant that led to the correlations 
of the form in eqn. (^). 

IV. SCALING AMPLITUDES 

The contribution of the short- distance terms may now be subtracted from the high- and 
low-temperature series, leaving the long-distance part, from which the scaling amplitudes 
may be computed using any of a variety of series analysis techniques. Such analysis is vastly 
simplified by the observation that there are no logarithmic or non-integer power contributions 
to the scaling-amplitude functions F±. 



To show this, independently of any fitting procedure, we have noted that any contribution 
to F± which is not a positive integer power of r would manifest itself in the high order 
series coefficients of the scaled susceptibility, (1 — s'^^)^^^^x±- The 1 + r/2 terms in F-t, as 
poles in the scaled susceptibility, also contribute but as their amplitudes are known to high 
precision, they can be subtracted. The residual coefficients are comparable in magnitude 
to those expected from the first neglected short-distance term which enters at r^^. We may 
place limits on the size of the amplitudes of any putative non-analytic terms in the scaling- 
amplitude functions. For example, for terms of the form ApT^ln |r|, the bounds, 

\Ap\ < 10-^^3007r(p-l), T>Tc, (7) 

\Ap\ < 10-^^6007r(p-l), T<T^, (8) 

reasonably exclude all p less than about 15. 

On purely numerical grounds, the absence of logarithmic corrections is surprising since it 
implies the cancellation of the many logarithmic multipliers in the scaling terms we discarded 
in the previous section. On the other hand, the absence of logarithms appears to be a 
requirement of the combination of the fermionic constraint p > q^ in (^ and renormalization- 
group scaling as discussed in the next section. 

To compute the amplitudes of the integer powers of r in the scaling-amplitude functions 
F±, we have carried out two independent analyses, one in the s-plane, the other in the 
f-plane, where v = tanhi^ is the conventional high-temperature variable. The natural 
boundary singularities at |s| = 1 are mapped to two circles, |f ±1| = v2- As they are farther 
from the origin than the ferromagnetic and antiferromagnetic singularities at f = ±(v2 — 1) 
their amplitudes are exponentially damped and may be neglected in the analysis. The s- 
plane analysis must take account of these singularities explicitly. The two analyses are in 
complete agreement. 

We find numerically that the scaling-amplitude functions multiphed by y/s appear to be 
even functions of r, the amplitudes of the odd terms being comparable in magnitude to 
the uncertainties in the even amplitudes. The rational amplitudes of r^ and r^ below we 
conjecture to be exact, and these values were fixed in the final fitting. The results, with 



uncertainty only in the final digits, are 

v^F+ = 1 + rV2 - rVl2 - 0.1235292285752086663r^ 
+ 0.136610949809095r^ - 0.13043897213r^° 
+ 0.1215129r^2 _ o.ll3ri^ + 0(r^^), 

V^F_ = 1 + rV2 - rVl2 

- 6.321306840495936623067r^ 

+ 6.25199747046024329r^ - 5.6896599756180r^° 

+ 5.142218271r^2 - 4.67472r^^ + 0(r^^). (9) 

V. COMPARISON WITH SCALING PREDICTIONS 



Prior to the analysis of [^ , all known amplitudes were in agreement with the hypothesis 
that corrections to scaling were due to scaling-field nonlinearity, and not to the presence of 
irrelevant variables. Here for the first time, we quantify the error in this "simple" scaling 
theory. Ignoring irrelevant operators, the expressions for the free energy, magnetization and 



susceptibility in zero magnetic field are [14 



/(r) = -A(ao(r))2ln|ao(r)|+Ao(r), 
M(r<0) = ml(r)|ao(r)|l/^ 
(3-\±{r) = C^{h{r)f\a,{T)\-'/^ 

-Ea2(r)ao(r) In |ao(r) I +Z}(r), (10) 

where A, B, C±, and E are constants and Ao{t) and D{t) are analytic functions of r. The 
functions ao(T) and 6i(t) are the leading terms in the expansion of the scaling fields, and 
can be determined from the free energy and magnetization. The result for x± is of the form 
(H) but with F replacing F± where 

JsF = l + + + 0r*. 11) 

^ 2 384 3072 ^ ^ ^ ^ ^ 

Note that this expression should hold in both temperature regimes. 

The difference between (|ll]) and @ we believe to be due to the effects of one or more 
irrelevant operators, confluent with the "simple" scaling terms. As there is no free param- 
eter to vary in the model, we can't identify these operators from the information we have. 

8 



However, it is likely that there are at least two mechanisms at work, one entering at O(r^) 
which preserves the equality of -F+ and F_, and a second entering at O(r^) which breaks 
this symmetry. In order to probe these effects further, we hope to study the model with 
anisotropy, and on other periodic lattices. 

The corrections to scaling we have found are confluent with expected analytic terms and 
in the renormalization group picture of scaling this leads to the possibility of logarithmic 
terms as well (cf. fl^). Logarithmic corrections are not demanded — the issue is whether 
the scaling fields are coupled and this depends on microscopic details. Barma and 



Fisher |I9| have investigated a model renormalization group flow in detail and conclude 
that in the case of a confluence, here labelled by integer index m, one must expect 
either no coupling between fields or corrections of the form (r^log |r|)'^ to all order 
k. Since the latter violates the fermionic constraint mk > k'^ we conclude there cannot 
be any logarithmic terms in the scaling-amplitude function F± as we have verified to 0(r^^). 
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